The stability of nonequilibrium polariton superflow in the presence of a cyhndrical 

defect 
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We make a theoretical study of the stabihty of nonequihbrium polariton superflows that interact 
with a cylindrical defect. The nonresonantly pumped polariton condensate is modelled with a 
generalized complex Ginzburg-Landau equation. At low pump intensities the dissipation is found 
stabilize the superflow. At large pump intensities, we find an instability that sets a lower critical 
speed for superfluidity. For even larger pump power, the lower and upper critical speed meet and 
stable superflows are no longer possible. 
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I. INTRODUCTION 

The superfluid properties of exciton polariton quan- 
tum fluids is a topic of active experimental and theoret- 
ical investigation. As predicted by Carusotto an Ciuti^, 
recent experiments by Amo et ali^ have shown that the 
scattering of the polaritons off a weak defect is strongly 
suppressed when the polariton speed is below the Landau 
critical speed. More recently, the scattering with a large 
defect was addressed theoreticallyS. and experimentally 
under continuous wave^ and pulsed^ excitation. In con- 
trast to the case of weak defects, where the supercurrents 
emit sound waves, strong defects lead to the nucleation 
of vortices and solitons. 

An important aspect of polariton condensates is the fi- 
nite polariton life time, that makes a profound difference 
with equilibrium Bose-Einstein condensates of dilute va- 
pors of ultracold atomic gases. When the finite polariton 
life time is compensated by continuous resonant excita- 
tion, the driving laser fixes the phase of the polariton 
fluid, which prevents the formation of phase defects such 
as vortices. In order to circumvent this phase rigidity, 
the defect was placed behind the excitation spot^. An 
alternative possibility to replenish the polariton conden- 
sate without fixing the phase of the condensate is by so- 
called nonresonant excitation. High energy excitations 
are injected by the laser and during their relaxation to- 
ward the bottom of the lower polariton branch, they loose 
all phase information and spontaneously form a coherent 
state when the density exceeds a threshold^. 

Where the superfiuid properties under resonant excita- 
tion are remarkably close to the ones of equilibrium con- 
densates, the nonresonant driving introduces profound 
differences. Because the decay rate is not negligible with 
respect to the interaction energy, the driven-dissipative 
nature of the polariton gas strongly affects the excitation 
spectrum^'* of polariton condensates and hence their su- 
perfluid properties. In a previous work, we have studied 
the faith of polariton superflows when they scatter with 
weak defects^. We have found that the dissipation actu- 
ally stabilizes the superflow and that the critical velocity 
is higher than in an equilibrium condensate at the same 
density. Even when excitations are created and a drag 



force is exerted on the defect, the condensate is not de- 
stroyed, but is able to maintain its superfluid velocity. 
A flnite concentration of defects is needed to destroy the 
superflow. 

From the theoretical side, much of the behavior of equi- 
librium superfluids can be understood from the Gross- 
Pitaevskii equation (GPE) that describes the order pa- 
rameter of a quantum degenerate Bose gasi^. In partic- 
ular, the GPE provides a description of the onset of a 
drag force on a defect that moves in a superfluid that 
sets on only when the defect defect velocity exceeds a 
critical speed Vc- For a weak defect, the GPE description 
reproduces the Landau criterion Vc — c, where c is the 
speed of sound-'^^. Above the critical speed, the emission 
of phonons is responisible for a drag force on the defect. 

The Gross-Pitaevskii equation is not limited to the 
limit of weak defects and allows for the investigation of 
nonperturbative defects that are dragged through a su- 
perfluid. Carrying out this analysis, Frisch et al. found 
that the critical velocity of a large cylinder is approx- 
imately two times smaller than the speed of sound 
In that case, the drag force on the defect results from 
the emission of vortex anti- vortex pairs^'^, that were re- 
cently observed in a Bose Einstein condensate of ultra- 
cold atomsi^. Above the critical speed, the continuous 
generation of excitations by the defect heats up the Bose 
gas, resulting in a final state of the gas moving along with 
the defect. 

The Gross-Pitaevskii equation can be straightfor- 
wardly generalized to nonequilibrium polariton conden- 
sates, where the steady state is determined by a bal- 
ance between driving and dissipation. Supplementing 
it with dissipative terms leads to the so-called complex 
Ginzburg-Landau equation (cGLE)i^, first used in the 
context of polariton condensation in Ref. 13 

In this work, we use the cGLE to analyze the colli- 
sion of a polariton condensate with a large diameter high 
cylindrical potential. For small dissipation, we find that 
the equilibrium picture continues to hold qualitatively: 
the superflow decays above a certain critical velocity. In 
analogy with the case of a perturbative defect, the critical 
velocity is increased by the dissipation. Surprisingly, we 
will find at that large dissipation, superfiuidity is com- 
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pletely destroyed by a strong defect. The reason is a non- 
trivial pattern of the currents in the condensate. In an 
intermediate window of dissipation, superfiows are stable 
between a lower and upper critical velocity. 

The outline of the paper is as follows. In Sec. HTJ we 
recapitulate the generalized complex Ginzburg Landau 
Equation (cGLE) that we use for the analysis of nonequi- 
libriuni superfluidity. In Sec. lIIIl we describe the conden- 
sate state in the presence of a large cylindrical defect. 
The stability of supeflows is then investigated in Sec.lIVI 
Conclusions are finally drawn in Sec. |Vl 



II. THEORETICAL MODEL 

We use in this work a generalization of the cGLE, that 
we have introduced in Ref. 1^. It includes a frequency- 
dependence of the pumping: 
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The energy zero has been set for convenience at the 
bottom of the lower polariton branch; the efficiency of 
amplification (proportional to the pumping strength P) 
decreases to zero in a frequency interval fl^ above it. 
Assuming a linear form of the frequency dependence of 
the amplification, the generalized GPE ([1]) maintains a 
temporally local form. The others terms describe gain 
saturation (r), losses (7), polariton mass (m), polariton- 
polariton interactions (g), external potential (Vext)- 

If the strength of pumping is large enough to overcome 
losses P > 7, the f/' = state is dynamically unstable 
against the creation of a finite condensate amplitude in 
any of the the low- momentum modes for which hk'^ /2m < 
VIk (1 — l/P)- The rate of this instability is maximum 
at = and decreases for increasing k. Under these 
conditions, a condensate is formed at fc = 0. 



III. CONDENSATE STATE 



When dissipation is included, this simple state can be 
dramatically modified. Steady states of the condensate 
state for two values of the pump power are shown in 
Fig. [TJ It is striking for the larger pump intensity, the 
condensate does no longer form close to zero momentum, 
but instead on a circle in momentum space. The arrows 
in the real space picture show the fiow of the condensate 
and indicate that the condensate fiows away from the 
defect. The physical explanation for the condensate flow 
is the following. The strong defect forces the density to 
vanish. Close to the defect, the density is below the value 
where the gain is saturated. Stabilization can then only 
occur when the excess created particles are carried away 
by a particle flow. This flow forces the condensation to 
occur at finite momentum even far away from the defect. 
This interpretation supports the numerical observation 
that the radius of the ring in momentum space increases 
with increasing pump power. 

The finite momentum states in Fig. [T] present flows of 
the condensate polaritons that are the gradient of the 
superfluid phase Vg — V9 and one may be tempted to 
conclude that the observation of such flows would prove 
superfluidity. The analogy of this type of superfiows 
with superfiuidity phenomena in equilibrium Bose gases 
is however difficult to sustain and strongly relies on the 
fact that energy is fed into the system. It is therefore 
analogous to the flow of water with a pressure difference 
between both ends. 





In the presence of an external potential Vexti the con- 
densation phenomenology becomes much richer. If Vext is 
a random potential, that naturally occurs because of the 
fluctuations during the microcavity growth, the forma- 
tion of multiple condensates was observed and explained 
within a related theoretical formalismii. 

In this work, we will study the effect of a cylindrical 
defect on the condensate. The amplitude is taken to be 
much larger than the polariton blue shift g|V'|^ due to 
the polariton-polariton interactions. In the absence of 
driving and dissipation, this type of potential makes a 
hole in the condensate that recovers to the asymptotic 
density on a distance fixed by the healing length ^ — 
^ h? /mgric, where ric is the condensate density. 



FIG. 1: Real (upper panels) and momentum (lower panels) 
space of a spontaneously formed polariton condensate in the 
presence of a defect. Left hand panels - Weakly out of equi- 
librium (small pumping and gain saturation), the condensate 
forms at A; = 0. Right hand panels ~ far from equilibrium, 
the condensate forms on a ring in momentum space, because 
of the currents away from the defect. Paramters: P/7 = 2 
(left) and P/7 = 15 (right); r/g = i, defect with diameter of 
and hight of IO7. Lengths and momenta are expressed in 
terms of — \Jh? j (1717) and fc^ — V^^ respectively. 
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FIG. 2: Critical velocity of a nonequilibrium condensate flow- 
ing against a big defect as a function of the gain saturation 
parameter r. Diamonds denote the standard critical velocity 
for breaking due to vortex emission, where the squares denote 
the lower critical velocity, below which the condensate state is 
unstable with respect to the formation of a ring in momentum 
space (see Fig. [3|. Parameters: g\ip'^ — 1.57, same defect as 
in Fig. m 



IV. STABILITY OF SUPERFLOW 

A true manifestation of superfluidity is the possibility 
of manifesting multiple flow patterns. In the absence of 
defects, flow of equilibrium superfluids is stable at any 
velocity, thanks to Galilean invariance. When a weak 
defect potential breaks this symmetry, the superflow loses 
stability above the speed of sound = ^Jgn^Jrn. For a 
large cylindrical defec t, the critical velocity is smaller and 
equalsi^ Wc = \/2/llc. 

In polariton condensates, condensation can be forced 
to occur in a momentum state different from zero by 
using a strong seed pulse at a finite wave vector. This 
acts as an initial condition for the condensate field and 
the stability of the superflow will determine whether the 
condensate continues to flow at that speed or will rather 
relax to a lower momentum. For simplicity, we work 
with periodic boundary conditions so that the flowing 
condensate does not leave the region where the cylinder 
is present. In practice this geometry could be realized 
with a ring shaped trapping potential. 

In our previous work^, we have investigated the stabil- 
ity of superflow for driven-dissipative condensates. We 
have pointed out that the energy dependence of the gain 
term breaks Galilean invariance and that therefore the 
critical velocity is finite even in the clean system. The 
effect of defects is on the other hand weakened by the dis- 
sipation. The critical velocity may be well above the one 
predicted for the equilibrium condensates Vc = \J gric/m. 

We will now turn to the investigation of the critical 
velocity when the nonequilibrium superfluid collides with 
a large cylinder. We have found that a single strong 
defect can lead to the nucleation of vortices and hence to 
the decay of the superflow. In order to assess the effect 



of the nonequilibrium situation on the stability of the 
superflow, we plot with diamonds in Fig|5] the critical 
velocity as a function of the parameter r / keeping the 
condensate density constant. Let us first concentrate on 
small values of the gain saturation parameter r/g. For 
r/g tending to zero, we recover the equilibrium critical 
velocity Vc ~ 0.5c. The critical velocity increases with 
the dissipation parameter r/g. In analogy to the case of 
weak defects, the dissipation stabilizes the superflow. 
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FIG. 3: Time evolved nonequilibrium condensate colliding 
with a big defect for various initial wave vectors, with in- 
creasing initial velocity from left to right. The upper panels 
show the real space density and flow (arrows), the lower panels 
show the momentum distribution. The squares indicate the 
initial wave vector of the condensate. Only when the initial 
velocity overcomes a critical value, the condensate maintains 
its momentum. For low initial velocity, a flow away from the 
defect appears. 

When the gain saturation is increased, the pump power 
should also increase to keep the density constant. At a 
certain value of the gain saturation (r/g w 1.5 for the 
density considered here), the condensate starts to form 
on a ring in momentum space with radius kc, as discussed 
above. From our numerical simulations, we have found 
that superflows are only stable if their wave vector is 
larger than kc- We thus find a lower critical velocity for 
the superflow, that is indicated with squares in Fig. [21 
The region of stable superflow is the shaded area. 

For large values of r/g, the upper critical velocity de- 
creases with gain saturation. At a critical value of r/g 
(w 3.5 for the considered density), the upper and lower 
critical velocity coincide and a stable superflow can no 
longer be generated. We assume that the physical reason 
is the following. The speed in the wake behind the defect 
is larger than the asymptotic speed that is imposed by 
the initial condition. For increasing gain saturation (and 
pumping so to keep the density fixed), the local velocity 
increases and therefore exceeds the critical velocity for 
lower values of the asymptotic velocity. 

This analysis has highlighted the role of gain satura- 
tion on the superfluid behavior of nonequilibrium conden- 
sates. In experiments however, it is rather the pump in- 
tensity than the gain saturation that is varied. In Fig. 31 
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FIG. 4: Lower (squares) and upper (diamonds) critical veloc- 
ity of a nonequilibrium condensate as a function of the pump- 
ing P. The gain saturation r = 0.3 and interaction strength 
<7 = 0.1. 



we show the lower and upper critical velocity as a func- 
tion of the pump power, for fixed r and g. It shows 
qualitatively the same trend as Fig. [2] The increase of 
the upper critical velocity with pump power is now com- 
ing from both the increased role of interaction energy 
gli/ip and gain saturation rjT/ip. The decrease of the up- 
per critical velocity at high powers is due to the effect 
of dissipation described above. The lower critical veloc- 
ity is zero for small powers {P < 1), increases and then 
saturates until it meets the upper critial speed. 

The enhanced stability of the superflow in nonequi- 
librium condensates with respect to the equilibrium ones 
raises the question whether the nonequilibrium superflow 
can be stable even without interactions between the po- 
laritons {g = 0). Fig. [5] shows that this is indeed the 
case. A stability region between the lower and upper 
critical velocity exists for pump intensities not too much 
above threshold. Gain saturation alone thus provides 
a nonlinearity that is sufficient to make superflows sta- 
ble agains scattering off defects. This observation is in 
stark contrast with equilibrium condensates where the 
only nonlinearity comes from interactions and therefore 
the noninteracting Bose gas has a zero critical velocity. 



with a cylindrical defect. The combination of the driven 
dissipative nature of polariton condensates with the de- 
fect potential can have a dramatic effect on the conden- 
sate shape. Above a critical pump power, the conden- 
sate is formed on a ring in momentum space, rather than 
around k = 0. When a superflow is imposed on the con- 
densate, a rich phenomenology is found. For small pump 
power/gain saturation, the equilibrium phenomelogy is 
recovered: stable superflow below a critical velocity and 




FIG. 5: The same as Fig.[4l but with zero interaction strength 
ff = 0. 



shedding of vortices above it. For increasing driving, the 
k = state becomes unstable and superflows are only 
stable above a lower critical velocity. At large driving, 
the upper critical velocity decreases and finally meets the 
lower one at the critical power above which no longer a 
stable superflow exists. Finally, we have shown that even 
in the absence of polariton-polariton interactions, a finite 
critical speed exists thanks to the gain saturation non- 
linearity. 
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